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Impedance  measurements  were  conducted  to  gain  insight  into  flooding  of  a  single  polymer  electrolyte 
membrane  (PEM)  fuel  cell.  The  stochastic  character  of  the  formation  of  water  droplets  and  subsequent 
removal  by  gas  flow  is  demonstrated  to  increase  the  standard  deviation  of  impedance  measurements, 
yielding  a  sensitive  manner  to  detect  onset  of  flooding.  The  increase  in  stochastic  noise  associated  with 
flooding  was  more  apparent  at  low  frequencies,  due  to  the  closer  match  to  the  characteristic  frequency 
associated  with  the  growth  and  removal  of  water  droplets.  The  onset  of  flooding  was  sensitive  to  the 
design  of  the  gas  diffusion  layer. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  performance  of  the  fuel  cell  is  influenced  by  kinetic  limita¬ 
tions  at  low  current  densities,  Ohmic  limitations  at  intermediate 
current  densities,  and  mass-transfer  limitations  at  high  current 
densities.  Kulikovs ky  [1  ]  and  Berg  et  al.  [2]  have  described  in  detail 
the  critical  role  of  water  management  in  operation  of  the  fuel  cell.  To 
maintain  proton  conductivity,  the  fuel  cell  membrane  must  remain 
hydrated.  To  achieve  hydration,  the  relative  humidity  of  inlet  gasses 
is  typically  held  at  a  large  value.  Water,  however,  is  also  a  product 
of  the  cathodic  reaction;  thus,  an  excess  of  water  in  the  cathode  is 
commonly  observed,  which  can  lead  to  condensation  and  subse¬ 
quent  flooding.  Flooding  increases  the  resistance  associated  with 
the  gas  diffusion  layer  and  may  even  block  flow  channels,  reducing 
the  availability  of  oxygen  [3].  Condensed  water  may  be  removed 
by  gas  flow.  Thus,  changes  in  design  of  reactant  flow  channels  and 
gas  diffusion  layers  have  been  proposed  to  reduce  the  influence  of 
flooding. 

A  comprehensive  review  of  the  methods  used  to  directly  mea¬ 
sure  the  presence  of  liquid  water  during  the  fuel  cell  operations 
was  provided  by  St-Pierre  [4].  Pressure  drop  has  been  reported 
to  provide  a  suitable  diagnostic  tool  for  monitoring  flooding  in 
the  fuel  cell  [5].  Flooding  was  also  investigated  by  correlating  the 
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appearance  of  flooding  to  the  Faradaic  resistance  [6].  Barbir  et  al. 
[7]  have  investigated  the  relationship  between  pressure  drop  and 
cell  resistance  to  make  a  distinction  between  flooding  and  drying. 
They  observed  that  both  pressure  drop  and  cell  resistance  changed 
in  case  of  drying;  whereas,  only  pressure  drop  changed  under 
flooding  conditions.  Ge  et  al.  [8]  observed  that  the  anode  flood¬ 
ing  is  mainly  due  to  water-droplet  condensation  at  channel  walls 
in  contrast  to  flooding  at  the  cathode  which  is  usually  attributed  to 
condensation  in  gas  diffusion  layer  (GDL).  They  have  also  reported 
that  use  of  a  hydrophilic  GDL  and  elevated  anode  plate  tempera¬ 
ture  could  mitigate  anode  flooding.  The  onset  of  flooding  may  be 
seen  in  steady-state  measurements,  but  the  impedance  response 
is  even  more  sensitive  to  appearance  of  flooding  conditions.  The 
impedance  technique  has  recently  been  used  to  detect  membrane 
drying,  flooding,  and  anode  poisoning  of  fuel  cell  stacks  [9].  Merida 
et  al.  [  10]  have  also  investigated  failure  modes  (drying  and  flooding) 
of  the  fuel  cell  using  the  impedance  technique.  The  approach  taken 
by  LeCanut  et  al.  [9]  and  Merida  et  al.  [10]  was  to  detect  flooding 
by  observing  increases  in  the  value  of  the  impedance  as  compared 
to  a  normal  impedance  measured  at  the  beginning  of  cell  opera¬ 
tion.  The  advantage  of  their  approach  is  that  a  physical  model  is 
not  needed.  The  difficulties  with  normalizing  the  impedance  to  the 
impedance  measured  at  the  beginning  of  cell  operation  are  that 
steady-state  operation  will  generally  not  be  established  during  this 
measurement,  that  there  may  be  other  reasons  for  increases  in  cell 
impedance  with  time,  and  that  flooding  may  already  be  taking  place 
during  the  initial  measurement. 
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Locally  resolved  impedance  spectroscopy  and  NMR  imaging 
have  been  used  to  investigate  flooding  and  drying  in  the  fuel  cell 
by  Schneider  et  al.  [11-13].  The  authors  reported  that,  for  co-flow 
configurations,  membrane  drying  was  evident  near  the  gas  inlet 
and  flooding  was  severe  near  the  gas  outlet  [12].  The  authors  have 
also  reported  that  drying  and  flooding  were  more  pronounced  in 
co-flow  as  compared  to  counter-flow  configurations  [13].  Fouquet 
et  al.  [14]  fitted  a  Randles-like  equivalent  circuit  to  impedance  data 
and  correlated  circuit  values  to  the  state-of-health  (flooding  and 
drying)  of  the  fuel  cell.  The  Randles-like  circuit,  however,  cannot 
account  for  all  the  phenomena  taking  place  in  the  fuel  cell.  While 
there  are  differences  in  the  specific  approaches  taken,  the  underly¬ 
ing  concept  for  each  of  these  approaches  was  that  one  can  detect 
flooding  by  observing  increases  in  the  value  of  the  impedance. 

The  object  of  the  present  work  was  to  explore  how  the  stochas¬ 
tic  character  of  flooding  can  be  exploited  to  improve  sensitive  of 
impedance  spectroscopy  to  detect  onset  of  flooding.  Impedance 
measurements  were  performed  as  a  function  of  different  parame¬ 
ters  such  as  current  density,  temperature,  back  pressure,  and  time. 
A  model  for  base-level  noise  in  impedance  measurements  for  nor¬ 
mal  conditions  (non-flooded)  was  developed  using  a  measurement 
model  analysis  [15-19],  and  stochastic  errors  were  also  assessed 
by  transient  fixed-frequency  measurements.  A  comparison  of  the 
actual  noise  to  the  base-level  noise  was  used  to  detect  onset  of 
flooding.  A  preliminary  version  of  this  concept  was  described  by 
Roy  and  Orazem  [20]. 

2.  Experimental 

The  experimental  system  and  the  impedance  instrumentation 
used  are  presented  in  this  section. 

2.1.  Materials  and  chemicals 

The  membrane  electrode  assembly  (MEA)  (purchased  from  Ion 
Power  Inc.,  New  Castle,  DE)  employed  0.0508  mm  (2  x  10-3in.) 
thick  Nation  N112  with  catalyst  layers  of  about  0.025  mm  on  both 
sides  of  the  membrane.  The  active  surface  area  was  5  cm2.  The  cat¬ 
alyst  layers  were  platinum  supported  on  carbon  with  a  Pt  catalyst 
loading  of  0.4  mg  cm-2  on  both  the  anode  and  the  cathode  sides. 
Two  types  of  gas  diffusion  layers  (GDL)  were  used  during  assem¬ 
bly  of  the  MEA.  Both  have  an  effective  thickness  of  0.284  mm  and 
were  made  of  carbon  cloth,  but  one  was  uniformly  macro-porous 
while  the  other  had  variable  porosity.  The  non-uniform  GDL  was 
micro-porous  to  the  catalyst  side  and  macro-porous  to  the  channel 
side.  Similar  GDL  structures  have  been  reported  in  the  literature 
[21-23].  The  material  of  the  interdigitated  flow  channel  used  was 
graphite  with  the  outlet  lower  than  the  inlet  to  facilitate  removal 
of  condensed  water.  A  torque  of  45  in.-pounds  was  applied  to  the 
fuel  cell  assembly.  Hydrogen  gas  was  used  as  fuel  and  a  79%  N2and 
21%  02 mixture  was  used  as  oxidant.  Compressed  N2  was  used  to 
purge  the  fuel  cell  before  and  after  experiments.  A  Barnstead  E-Pure 
Water  System  with  an  ion  resistivity  of  14.9  cm  was  used  as  a 
source  of  deionized  water  delivered  to  the  anode  and  the  cathode 
humidifiers. 

An  850C  fuel-cell  test  station  (supplied  by  Scribner  Associates, 
Southern  Pines,  NC)  was  used  to  control  reactant  flowrates  and 
temperatures.  The  test  station  was  connected  to  a  computer  by 
an  interface  for  data  acquisition.  The  gas  flows  were  humidified  to 
100%  relative  humidity  at  the  respective  temperatures.  The  hydro¬ 
gen  flow  rate  was  0.1 1  min-1  and  the  air  flow  rate  was  0.5 1  min-1 . 
The  maximum  stoichiometry  for  hydrogen  and  air  was  1.5  and  2.5, 
respectively  and  the  cell  was  operated  at  the  fully  humidified  con¬ 
dition. 


Fig.  1.  Impedance  data  recorded  with  the  850C  with  applied  current  density  as  a 
parameter.  The  anode,  cathode,  and  cell  temperatures  were  set  to  50  °C. 

2.2.  Electrochemical  impedance  measurements 

Impedance  measurements  were  performed  using  two  differ¬ 
ent  systems.  The  850C  fuel-cell  test  station  contains  an  electronic 
load  and  frequency  response  analyzer.  All  electrochemical  mea¬ 
surements  were  performed  with  a  two-electrode  cell.  The  anode 
was  used  as  a  pseudo-reference  electrode.  The  impedance  mea¬ 
surements  were  conducted  in  galvanostatic  mode  for  frequency 
range  of  10  kHz  to  5  mHz  with  a  10  mA  peak-to-peak  sinusoidal 
perturbation.  Lissajous  plots  were  used  to  confirm  that  the  pertur¬ 
bation  amplitude  was  low  enough  to  achieve  a  linear  response.  The 
corresponding  potential  perturbation  ranged  from  0.04  to  0.4  mV. 
The  frequencies  were  spaced  in  logarithmic  progression  with  10 
points  per  frequency  decade.  Impedance  scans  were  conducted  in 
auto-integration  mode  with  a  minimum  of  2  cycles  per  frequency 
measured. 

3.  Results 

The  influence  of  flooding  on  the  operation  of  the  fuel  cell  was 
investigated  using  impedance  spectroscopy.  The  results  of  fre¬ 
quency  scan  and  single-frequency  time-dependent  measurements 
are  presented  in  following  sections. 

3.1.  Impedance  response 

A  typical  impedance  response  is  presented  in  Fig.  1  with  current 
density  as  a  parameter.  The  size  of  the  intermediate-frequency  and 
the  low-frequency  arcs  increased  with  increasing  current  density, 
an  effect  which  was  suggested  in  the  literature  [9,14]  to  be  due  par¬ 
tially  to  flooding.  The  impedance  spectra  were  relatively  smooth  for 
low-current  densities;  however,  the  spectra  have  significant  scat¬ 
ter  at  higher  current  densities  where  flooding  was  probable.  The 
scatter  was  particularly  evident  at  low  frequency,  where  the  spec¬ 
tra  show  jumps  in  value  which  may  be  associated  with  removal  of 
condensed  water.  Similar  results  were  observed  at  other  tempera¬ 
tures.  The  impedance  measured  at  70 0  C  is  presented  in  Fig.  2  with 


Fig.  2.  Impedance  data  recorded  with  the  850C  with  applied  current  density  as  a 
parameter.  The  anode,  cathode,  and  cell  temperatures  were  set  to  70  °C. 
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Fig.  3.  Single-frequency  impedance  measurements  recorded  at  0.1  Hz,  70  °C,  and 
1.4  A  cm-2  as  functions  of  time:  (a)  real  part;  (b)  imaginary  part. 

current  density  as  a  parameter.  Here  also  the  impedance  data  have 
greater  degree  of  scatter  at  higher  current  densities,  and  this  scatter 
is  particularly  evident  at  lower  frequencies. 

3.2.  Stochastic  error  in  impedance  response 

The  standard  deviations  of  the  stochastic  errors  were  estimated 
from  both  impedance  spectra  and  single-frequency  transients  for 
dry,  flooded,  and  non-flooded  conditions. 

3. 2 A.  Sensitivity  to  flooding 

The  impedance  was  recorded  as  a  function  of  time  for  different 
current  densities  and  fixed-frequencies.  For  example,  the  real  and 
imaginary  parts  of  impedance  at  0.1  Hz  and  1.4  A  cm-2  are  pre¬ 
sented  as  functions  of  time  in  Fig.  3  (a)  and  (b),  respectively.  The 
standard  deviations  in  the  impedance  data  were  calculated  using 
a  moving  average  method  to  account  for  the  systematic  changes 
shown  in  Fig.  3.  The  standard  deviations  for  the  real  and  imagi¬ 
nary  parts  of  the  impedance  response  are  presented  in  Fig.  4(a)  for 
a  measured  frequency  of  100  Hz  and  in  Fig.  4(b)  for  a  frequency 
of  1  Hz.  At  low  frequencies,  Fig.  4(b),  the  standard  deviation  of  the 
real  part  of  the  impedance  is  clearly  larger  than  that  of  the  imagi¬ 
nary  part.  The  solid  line  given  in  Fig.  4  represents  the  model  value 
for  the  standard  deviation,  developed  in  a  subsequent  section  for 
non-flooded  conditions. 

The  standard  deviations  for  the  real  part  of  the  impedance  are 
presented  as  a  function  of  current  density  in  Fig.  5  for  frequencies  of 
0.1, 1,  and  100  Hz.  The  stochastic  errors  in  the  impedance  increased 


Fig.  4.  The  standard  deviation  of  the  single-frequency  impedance  measurements 
recorded  at  1.4  A  cm-2  and  70  °  C  as  functions  of  time:  (a)  at  a  frequency  of  100  Hz;  (b) 
at  a  frequency  of  1  Hz  (as  presented  in  Fig.  3).  The  solid  line  represents  the  empirical 
model  developed  for  the  error  structure  given  by  Eq.  (1). 

with  increasing  current  density,  in  particular  for  the  lower  frequen¬ 
cies  of  0.1  and  1  Hz.  The  standard  deviation  at  lower  frequencies  are 
higher  than  that  at  higher  frequencies  for  any  given  current  density. 

The  statistical  nature  of  the  error  structure  may  be  used  to 
explain  the  larger  standard  deviation  observed  at  larger  current 
densities.  As  shown  in  Fig.  4(b),  the  standard  deviation  of  the  real 
part  of  the  impedance  at  1  Hz  was  much  higher  as  compared  to 


Fig.  5.  The  standard  deviations  for  the  real  part  of  the  impedance  as  a  function  of 
current  density  with  frequency  as  a  parameter  for  cell  operation  at  70  °C. 
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Fig.  6.  Standard  deviations  for  the  impedance  data  obtained  at  a  current  density 
of  0.4  A  cm-2.  The  solid  line  represents  the  empirical  model  developed  for  the  error 
structure  given  by  Eq.  (1 ).  The  dashed  lines  represent  the  asymptotic  behavior  of  the 
model  at  high  and  low  frequencies. 


imaginary  part.  The  results  presented  in  Fig.  4(a)  indicate  that 
at  100  Hz,  the  standard  deviation  of  real  and  imaginary  parts  of 
the  impedance  were  also  unequal,  but  the  imaginary  part  of  the 
impedance  seems  to  have  the  larger  standard  deviation.  For  causal 
systems  which  satisfy  the  Kramers-Kronig  relations,  the  standard 
deviation  in  the  real  and  imaginary  parts  should  be  equal  [24].  The 
literature  indicates  that,  at  low  frequencies,  the  real  part  of  the 
impedance  is  more  sensitive  to  flooding  [9,10] ;  therefore,  the  higher 
standard  deviation  of  the  real  part  of  the  impedance  observed  at  low 
frequencies  may  be  attributed  to  onset  of  flooding. 


3.2.2.  Baseline  error  structure 

To  establish  a  baseline  error-structure  model  for  the  standard 
deviation  of  impedance  measurements  in  the  absence  of  flooding, 
a  measurement  model  analysis  [15-17,19]  was  applied  to  a  large  set 
of  replicated  impedance  data.  The  measurement  model  was  used 
to  filter  small  systematic  changes  from  one  measurement  to  the 
other.  The  standard  deviations  for  impedance  response  recorded  at 
0.4  A  cm-2,  presented  in  Fig.  6,  were  smaller  than  those  observed 
at  both  larger  and  at  smaller  current  densities.  The  real  and  imagi¬ 
nary  parts  of  the  impedance  were  statistically  indistinguishable  at 
all  frequencies,  in  agreement  with  expectations  for  data  that  are 
consistent  with  the  Kramers-Kronig  relations. 

The  model  for  error  structure  developed  in  previous  work  for  a 
broad  variety  of  electrochemical  and  electronic  systems  [25,26]  did 
not  provide  a  good  representation  of  the  error  structure  presented 
in  Fig.  6.  The  failure  of  the  general  model  to  apply  to  the  present  data 
was  attributed  to  differences  in  the  parameters  used  to  make  the 
measurement.  An  empirical  model  was  found  to  follow  the  form: 

&r  =  Oj  =  C  +  \Zr^max\{of  k)  (1) 

where  a  =  9  x  10-5,  b  =  0.695,  and  c  =  3.5  x  10-5.  The  model  can 
be  rationalized  by  examination  of  the  standard  formula  for  prop¬ 
agation  of  stochastic  errors,  which  can  be  written  for  impedance 
as 


"  dz  " 

2 

7 

dz  " 

9  A  / 

+  °v 

9AV 

where  A /  and  AV  are  perturbations  in  current  and  potential, 
respectively,  and  Z  is  the  impedance  given  as 


Z  = 


AV 

~A[ 


(3) 


Fig.  7.  Impedance  measurement  recorded  at  0.4  A  cm-2 and  1  Hz  as  functions  of  time 
at  70  °C:  (a)  real  part;  (b)  imaginary  part. 
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For  galvanostatic  modulation  at  fixed  amplitude,  Eq.  (5)  takes 
the  form: 


oz  =  B\Z\  +  C  (6) 

where  C  and  B  are  constants. 

The  error  structure  model  described  by  Eq.  (1)  is  empirical 
and  provided  a  better  representation  of  the  experimental  error 
structure  than  did  Eq.  (6).  The  error  structure  is  dependent  on 
instrumental  parameter  settings  which  are  not  accounted  for  in 
the  development  of  Eq.  (6).  While  the  above  development  did  not 
yield  the  best  model  for  the  error  structure,  it  serves  to  show  that  a 
dependence  greater  than  first  order  on  the  impedance  value  is  not 
expected. 

The  error  structure  identified  for  impedance  spectra  can  be 
compared  to  that  obtained  for  single-frequency  transient  measure¬ 
ments  such  as  presented  in  Fig.  7  for  impedance  data  collected  at 
0.4  A  cm-2  and  1  Hz.  The  standard  deviations  for  the  real  and  imag¬ 
inary  parts  of  the  impedance  response  are  presented  in  Fig.  8(a) 
for  a  measured  frequency  of  100  Hz  and  in  Fig.  8(b)  for  a  fre¬ 
quency  of  1  Hz.  The  error  structure  model  represented  by  Eq.  (1) 
is  in  good  agreement  with  the  standard  deviations  obtained  by 
transient  single-frequency  measurements.  At  100  Hz,  Fig.  8(a),  the 
standard  deviation  of  the  real  and  imaginary  parts  of  the  impedance 
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Fig.  8.  The  standard  deviation  of  the  single-frequency  impedance  measurements 
recorded  at  0.4  A  cm_2and  70  °  C  as  functions  of  time:  (a)  at  a  frequency  of  100  Hz;  (b) 
at  a  frequency  of  1  Hz  (as  presented  in  Fig.  7).  The  solid  line  represents  the  empirical 
model  developed  for  the  error  structure  given  by  Eq.  (1 ). 


Fig.  9.  Standard  deviations  for  the  impedance  data  obtained  at  a  current  den¬ 
sity  of  0.4 Acm-2:  (a)  with  system  temperature  as  a  parameter;  (b)  at  70°  C  with 
anode/cathode  back-pressure  as  a  parameter.  The  solid  line  represents  the  empirical 
model  for  the  error  structure  given  by  Eq.  (1 ). 


are  equal;  however,  as  shown  in  Fig.  8(b),  the  standard  deviations 
for  real  and  imaginary  parts  are  not  equal  at  a  frequency  of  1  Hz. 
This  result  suggests  that,  even  at  a  current  density  of  0.4  A  cm-2, 
where  the  error  structure  was  the  smallest,  some  flooding  may  be 
taking  place. 

Eq.  (1)  provided  a  good  representation  of  the  impedance  error 
structure  obtained  at  a  current  density  of  0.4  A  cm-2  under  a 
broad  variety  of  conditions.  The  error  structure  for  the  impedance 
response  collected  at  different  temperatures  is  presented  in  Fig.  9(a) 
and  the  influence  of  anode/cathode  back-pressure  is  explored  in 
Fig.  9(b).  The  model  provided  a  good  representation  of  the  error 
structure  for  all  cases  considered  in  Fig.  9.  Accordingly,  Eq.  (1 )  was 
used  to  represent  the  base-level  standard  deviation  for  measure¬ 
ments  unaffected  by  drying  or  flooding  conditions. 

3.2.3.  Detection  of  flooded  conditions 

Impedance  spectra  were  obtained  at  different  current  densities. 
The  standard  deviations  obtained  at  larger  current  densities  are 
compared  in  Fig.  10  to  the  values  obtained  at  0.4  Acm-2.  The  stan¬ 
dard  deviations  obtained  for  current  densities  of  1.0  and  1.4  Acm-2 
are  clearly  larger  than  that  predicted  by  Eq.  (1).  The  discrepancy 
may  be  attributed  to  stochastic  processes  within  the  cell  such  as 
associated  with  flooding. 

The  standard  deviations  presented  in  Fig.  5  were  normalized  by 
the  base-stochastic  errors  level  calculated  using  the  error  structure 
presented  by  Eq.  (1).  The  result  for  the  real  part  of  the  impedance 
is  presented  in  Fig.  11.  As  shown  in  Fig.  11,  the  standard  deviation 


of  the  impedance  data  increased  with  an  increase  in  the  operating 
current.  The  increased  noise  levels  are  seen  at  a  frequency  of  100  Hz 
as  well  as  at  lower  frequencies.  The  normalized  standard  deviations 
for  low  current  densities  have  some  values  less  than  unity,  which 
may  be  attributed  to  uncertainty  in  the  model  for  the  baseline 
error  structure  due  to  the  inability  to  avoid  flooding  completely  in 
the  data  used  to  assess  the  baseline  error  structure.  At  large  current 
densities  and  at  low  frequencies,  the  standard  deviation  calculated 


f  /  Hz 

Fig.  10.  Standard  deviations  for  the  impedance  data  obtained  at  a  current  densities 
of  0.4, 1.0,  and  1.4  A  cm-2.  The  solid  line  represents  the  empirical  model  for  the  error 
structure  given  by  Eq.  (1 ). 
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Fig.  11.  Normalized  standard  deviations  for  the  real  part  of  the  impedance  calculated 
from  the  data  shown  in  Fig.  3  as  a  function  of  current  density  with  frequency  as  a 
parameter. 

for  the  real  part  of  the  impedance  was  more  than  20  times  that 
obtained  in  the  absence  of  flooding.  The  standard  deviation  of 
the  imaginary  part  of  the  impedance  had  no  clear  dependency  on 
flooding.  The  standard  deviation  of  the  real  part  particularly  at  low 
frequency,  however,  can  be  used  to  detect  onset  of  flooding. 

Similar  experiments  and  analysis  were  performed  on  the  MEA 
with  a  uniform  GDL.  The  impedance  response  is  presented  in  Fig.  12 
with  current  density  as  a  parameter.  The  scattering  at  higher  cur¬ 
rent  densities  was  more  evident  than  observed  for  experiments 
using  the  non-uniform  GDL.  Transient  single-frequency  impedance 
measurements  were  used  to  obtain  the  standard  deviations  for  real 
and  imaginary  parts  of  the  impedance.  The  standard  deviation  for 
the  real  part  of  the  impedance  is  presented  in  Fig.  13  as  a  function  of 
current  density  for  frequencies  of  0.1, 10,  and  100  Hz.  The  baseline 
error  structure  model  used  to  normalize  the  data  was  calculated  at 
a  current  density  of  0.2  A  cm-2. 

3.2.4.  Detection  of  dry  conditions 

The  stochastic  errors  for  small  current  densities  are  shown  in 
Fig.  14  for  the  MEA  with  a  non-uniform  pore  distribution.  The 
empirical  model  given  by  Eq.  (1)  provided  a  good  description  for 
the  behavior  at  a  current  density  of  0.4  A  cm-2,  but  the  observed 
errors  are  much  larger  for  lower  current  densities  where  dry  condi¬ 
tions  are  anticipated.  This  result  is  consistent  with  the  experimental 
observations  of  Schneider  et  al.  [12]  who  report  impedance  scans 
with  large  scatter  at  low  frequencies  for  dry  conditions.  Scatter 
observed  for  poorly  humidified  gas  flows  has  been  attributed  in  the 
literature  to  the  appearance  of  multiple  steady  states  [27,28].  In  the 
present  case,  however,  properly  calibrated  humidification  bottles 
were  to  used  to  humidify  the  gas  flows.  The  drying  at  low  current 


Fig.  12.  The  impedance  data  recorded  using  the  MEA  with  a  uniform  GDL.  The  anode, 
the  cathode,  and  cell  temperatures  were  set  at  50  °C. 


Fig.  13.  Normalized  standard  deviations  for  the  real  part  of  the  impedance  as  a  func¬ 
tion  of  current  density  with  frequency  as  a  parameter  for  the  MEA  with  a  uniform 
pore  distribution.  The  anode,  the  cathode,  and  cell  temperatures  were  set  at  50  °C. 

densities  can  be  attributed  to  reduced  production  of  water  at  the 
cathode  coupled  with  redistribution  by  electro-osmosis.  This  work 
shows  that  the  standard  deviation  of  the  real  part  particularly  at 
low  frequency,  however,  can  be  used  as  well  to  detect  the  presence 
of  dry  conditions. 

4.  Discussion 

Due  to  its  influence  on  mass  transfer  and  kinetics,  the  onset 
of  flooding  in  the  fuel  cell  can  be  identified  by  a  decrease  in  cell 
potential  at  fixed  current  or  a  decrease  in  current  at  fixed  potential. 
Similarly,  the  onset  of  flooding  can  be  identified  by  an  increase  in 
the  cell  impedance.  In  the  present  work,  the  increase  in  the  low- 
frequency  cell  impedance  associated  with  flooding  was  on  the  order 
of  20%.  The  increase  of  20%  in  the  cell  impedance  was  accompa¬ 
nied  by  an  increase  in  the  standard  deviation  of  the  real  part  of  the 
impedance  by  2000  or  3000%. 

The  increase  in  stochastic  errors  in  the  impedance  measurement 
can  be  attributed  to  the  random  character  of  the  flooding  process  in 
which  droplets  of  water  are  formed  and  then  removed  by  gas  flow. 
In  fact,  the  increase  in  stochastic  errors  provides  verification  that 
the  increased  cell  impedance  was  at  least  partially  due  to  flooding. 
Impedance  spectroscopy  has  been  shown  to  provide  a  more  sensi- 


f  /  Hz 

Fig.  14.  Standard  deviations  for  the  impedance  data  obtained  at  a  current  densities 
of  0.02,  0.1,  and  0.4  A  cm-2.  The  solid  line  represents  the  empirical  model  for  the 
error  structure  given  by  Eq.  (1 ).  The  anode,  the  cathode,  and  cell  temperatures  were 
set  at  70  °C. 


218 


S.K.  Roy,  M.E.  Orazem  /  Journal  of  Power  Sources  184  (2008)  212-219 


40  (“i — i - 1  i  |  i  |  i  |  i  | - 1 — | — i  |  i 

_  -  -A-  -  MEA  with  non-uniform  GDL 

35  - 

-  -o-  -  MEA  with  uniform  GDL 


30  -  ® 


0.0  0.2  0.4  0.6  0.8  1.0  1.2  1.4 

Current  Density  /  A  cm  2 

Fig.  15.  Normalized  standard  deviations  for  the  real  part  of  the  impedance  measured 
at  0.1  Hz  for  fuel  cells  containing  two  different  MEAs  as  a  function  of  current  density. 
The  experiments  for  the  uniform  MEA  were  performed  at  50  °C,  and  the  experiments 
for  the  non-uniform  MEA  were  performed  at  70  °C. 

tive  assessment  of  cell  condition  than  steady-state  measurements 
of  cell  potential  and  current.  The  difficulty  with  using  impedance 
directly  to  detect  flooding  is  that  a  baseline  value  for  the  impedance 
must  be  established  in  the  absence  of  flooding.  This  baseline  must, 
however,  change  with  time  due  to  systematic  changes  to  catalyst 
and  membrane  properties  that  are  not  associated  with  flooding. 
Thus,  a  baseline  established  when  a  cell  is  first  commissioned  will 
not  be  valid  throughout  the  lifetime  of  the  cell. 

In  contrast,  the  model  given  by  Eq.  (1)  for  the  standard  devia¬ 
tion  of  the  impedance  should  be  affected  largely  by  instrumental 
settings,  and,  so  long  as  the  impedance  is  measured  in  the  same 
way,  a  baseline  established  for  the  standard  deviation  of  impedance 
measurements  in  the  flooded  condition  should  be  valid  throughout 
the  lifetime  of  the  cell.  In  addition,  the  change  in  the  standard  devi¬ 
ation  of  the  measurement  caused  by  flooding  is  100  times  larger 
than  the  corresponding  change  in  the  value  of  the  impedance.  Thus, 
assessment  of  the  standard  deviation  of  impedance  measurements 
will  provide  a  more  sensitive  indicator  for  the  onset  of  flooding  in 
a  PEM  fuel  cell. 

The  standard  deviation  of  the  impedance  measurements  were 
sensitive  to  the  properties  of  the  MEA  used  in  the  experiment.  A 
comparison  of  the  normalized  standard  deviations  for  the  real  part 
of  the  impedance  is  presented  in  Fig.  15  with  GDL  properties  as  a 
parameter.  The  experiments  for  the  uniform  MEA  were  performed 
at  50  °C,  and  the  experiments  for  the  non-uniform  MEA  were  per¬ 
formed  at  70  °C.  The  normalized  standard  deviation  increased  at 
lower  current  densities  for  the  GDL  with  a  uniform  pore  distribu¬ 
tion.  The  value  was  closer  to  unity  over  a  broader  range  for  the 
GDL  with  a  non-uniform  pore  distribution.  The  results  are  consis¬ 
tent  with  the  observation  that  a  larger  maximum  current  density 
could  be  obtained  with  the  non-uniform  GDL.  Micro-macro-porous 
GDLs  are  reported  to  provide  better  water  management  [21-23]. 
The  increase  in  normalized  standard  deviation  at  low  current  den¬ 
sities  observed  for  the  non-uniform  GDL  is  likely  due  to  drying  of 
the  membrane  [7,9].  Similar  experiments  were  not  performed  for 
the  uniform  GDL. 

5.  Conclusions 

The  flooding  of  gas  diffusion  layer  pores  in  the  fuel  cell  has  been 
associated  with  increases  in  the  internal  cell  resistance  and  in  the 
impedance  response  of  the  fuel  cell.  The  formation  and  removal  of 
water  droplets  is  an  inherently  stochastic  process  which  increases 
the  stochastic  errors  observed  in  impedance  measurements.  A  mea¬ 


surement  technique  oriented  towards  assessment  of  the  stochastic 
errors  can  therefore  be  used  to  identify  the  onset  of  flooding.  In 
the  present  work,  impedance  spectroscopy  was  coupled  with  a 
measurement-model-based  error  analysis  to  detect  onset  of  flood¬ 
ing.  This  method  is  particularly  attractive  because  it  is  extremely 
sensitive  and  a  well-defined  baseline  stochastic  error  can  be  estab¬ 
lished  for  the  non-flooding  condition. 

The  onset  of  flooding  was  examined  for  a  5  cm2PEM  fuel  cell 
with  an  interdigitated  flow  channel.  At  low-current  densities,  the 
ratio  of  the  observed  standard  deviation  to  the  expected  non- 
flooded  standard  deviation  was  close  to  unity.  At  larger  current 
densities,  the  ratio  for  the  real  part  of  the  impedance  became  quite 
large,  with  onset  of  flooding  evident  at  current  densities  above 
1 A  cm-2  for  the  MEA  with  a  non-uniform  GDL  and  at  current  den¬ 
sities  above  0.3  A  cm-2  for  the  MEA  with  a  uniform  GDL.  Drying  was 
evident  at  current  densities  below  0.3  A  cm-2  for  the  MEA  with  a 
non-uniform  GDL. 

The  work  presented  here  demonstrates  that  the  stochastic  error 
structure  of  impedance  measurements  may  be  used  to  detect  oper¬ 
ating  conditions  of  the  fuel  cell  which  induce  flooding  or  drying.  In 
this  case,  the  flooding  or  drying  phenomena  contribute  stochastic 
errors  which  are  superposed  on  those  associated  with  the  elec¬ 
tronic  instrumentation.  As  shown  in  the  present  work,  onset  of 
flooding  is  gradual  and  can  be  affected  by  factors  such  as  flow  chan¬ 
nel  design  and  GDL  characteristics.  The  current  density  reported 
in  the  present  work  for  onset  of  flooding  may  be  useful  for  the 
systems  investigated  here,  but  the  approach  is  general  and  can  be 
exploited  to  evaluate  flooding  and  drying  characteristics  for  larger 
fuel  cells  and  for  other  cell  designs.  The  next  step  in  the  develop¬ 
ment  of  this  technique  would  be  to  couple  the  impedance  error 
analysis  approach  with  methods  for  direct  measurement  of  liquid 
water  in  the  fuel  cell  such  as  neutron  imaging  [12]  or  the  residence 
time  distribution  technique  [29,30]. 
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